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Aluminum is the primary structural material in nuclear engineering, and its cross-section induced by 14 
MeV neutrons is of great significance. To address the issue of insufficient accuracy for the 7’Al(n,2n)7°Al 
reaction cross-section, the activation method and accelerator mass spectrometry (AMS) technique were used to 
determine the 7’Al(n,2n)*°Al cross-section, which could be used as a D-T plasma ion temperature monitor in 
fusion reactors. At the China Academy of Engineering Physics (CAEP), neutron activation was performed using 
a K-400 neutron generator produced by the T(d,n)*He reaction. The *°Al/*’Al isotope ratios were measured 
using the newly installed GYIG 1 MV AMS at the Institute of Geochemistry, Chinese Academy of Sciences. 
The neutron flux was monitored by measuring the activity of °”"Nb produced by the *Nb(n,2n)””Nb reaction. 
The measured results were compared with available data in the experimental nuclear reaction database, and 
the measured values showed a reasonable degree of consistency with partially available literature data. The 
newly acquired cross-sectional data at 12 neutron energy points through systematic measurements clarified the 
divergence, which has two different growth trends from the existing experimental values. The obtained results 
are also compared with the corresponding evaluated database, and the newly calculated excitation functions with 
TALYS-1.95 and EMPIRE-3.2 codes, the agreement with CENDL-3.2, TENDL-2021 and EMPIRE-3.2 results 
are generally acceptable. A substantial improvement in the knowledge of the 7’Al(n,2n)”°Al reaction excitation 
function was obtained in the present work, which will lay the foundation for the diagnosis of the fusion ion 


temperature, testing of the nuclear physics model, and evaluation of nuclear data, etc. 
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I. INTRODUCTION 


Its high strength-to-weight ratio and relatively low activa- 
tion rate make aluminum one of the most attractive materi- 
als for many structural applications, such as fusion reactors, 
accelerator subcritical systems, and Generation IV reactors 
[1, 2]. The cross-section of aluminum is of great significance 
not only for its general applicability in testing nuclear models 
but also for the estimation of displacement damage in A1- 
metal alloys and the radioactivity of the long-lived ground 
state *°Al (T1 /2 =7.17 x 10° a) that produced 1808.65 keV 
gamma ray [3]. In addition, a high threshold (13.5 MeV) 
and a dramatic increase in the 27 Al(n,2n)°Al reaction cross- 
section with neutron energy in the 14 MeV region plays an 
essential role in the monitoring of ion temperature in a D- 
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T fusion plasma, as pointed out by Smither and Wallner et 
al. [4, 5]. 


However, owing to the long half-life of ?6A1(Ma) and the 
small cross-section (mb) of the 21 1(n,2n)°Al reaction near 
the threshold, it is challenging to quantify the cross-section 
of the 7’Al(n,2n)*°Al reaction using gamma spectrometry. 
As the most sensitive measurement technology for 26A] to 
date [6], accelerator mass spectrometers (AMS) have signifi- 
cantly shortened the irradiation and measurement time of the 
°741(n,2n)*°Al reaction channel and effectively improved the 
accuracy of the data. In particular, accelerator mass spec- 
trometers have been widely used to measure the cross sections 
and half-lives of nuclear reactions associated with long-lived 
nuclides [7-11]. Nevertheless, cross-sectional measurements 
of the 27 Al(n,2n) A] reaction are insufficient, and the avail- 
able data differ significantly, as illustrated in Fig. 1, because 
of the limitations of single-energy neutron sources and high- 
precision accelerator mass spectrometers. 


The following is a brief overview of the previous 
°7A1(n,2n)*°Al reaction measurement procedure. With the 
emergence of accelerator mass spectrometry technology, the 
first measurement of the 27 A1(n,2n) 26A] cross-section was re- 
ported by Smither et al. [4] in 1984; however, these data were 
not included in the experimental nuclear reaction data (EX- 
FOR). Using a well-shielded germanium detector, Sasao de- 
tected 1.809 MeV y rays following the 8 decay of 26A] to cal- 
culate the 21 Al(n,2n)”®Al reaction cross-section in 1987 [12]. 
Although the uncertainty of the neutron energy was mea- 
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Fig. 1. (Color online) Previous results for the excitation function of 
the 21 AlI(n,2n)”% A1 reaction near threshold. 


sured by the ratio of the 7 r(n,2n) 89m+87r cross-section to the 
°2Nb(n,2n)?””Nb near 14 MeV, the cross-section error details 
were not given. In 1988, Iwasaki et al. conducted irradiation 
for several months using a DT generator. The measurements 
were performed after approximately one year of cooling. In 
1991, 15-34 MeV neutrons were produced by bombarding 
a 1 mm-thick beryllium target with a 3-5 microamp proton 
beam. Nakamura et al. [13] provided only high-energy cross- 
sectional data using the AMS. In 1996, Al samples were irra- 
diated for four days, and y-ray measurements lasted for more 
than 1.7 years using a Ge detector with high detection effi- 
ciency. The uncertainty of the neutron energy is not given 
in the original literature, whereas the cross-sectional uncer- 
tainty is in the range of 10% to 50%. In 1998, two samples 
were irradiated with a D-T neutron beam for approximately 
20 d. Qiang et al. obtained the first °7A1(n,2n)*Al cross- 
section at two neutron energy points using a 6 MV AMS [14] 
in China, and the data were the maximum cross-section in 
the 14 MeV neutron energy region. In 2000, Sudbrock et 
al. reported four different reaction cross-section measure- 
ments at 14.6 MeV neutron energy via AMS [1]. Although 
the error of the 21 A1(n,2n)Al reaction in Sudbrock’s work is 
more than 70%, the cross-section also falls within the range 
of values predicted from the systematics of the (n,2n) re- 
action cross sections at 14.6 MeV. In 2003, Wallner et al. 
[5] provided a comprehensive study on the measurement of 
the 27 Al(n,2n)*°Al reaction cross-section. At 14.8 MeV, the 
cross-sectional values of the 27 A1(n,2n)%Al reaction are 30% 
and 40% lower than Y. Ikeda’s and Q. Zhao’s, respectively. In 
2016, A.A. Filatenkov gave results consistent with Wallner’s 
study. 


Although the °7A(n,2n)*°Al reaction was carefully mea- 
sured and evaluated, as shown in Fig. 1, there are still non- 
negligible discrepancies between the values recommended by 
the different evaluated nuclear data libraries. The 14 MeV 
neutron energy region has a 50% difference between the eval- 
uated values [15-18]. The evaluated values of CENDL-3.2 
were significantly different from those of CENDL-3.1. Con- 


sidering the insufficient experimental data, further investiga- 
tion is required for the evaluated database based on various 
theoretical models. In conclusion, there are two distribu- 
tion trends in the existing experimental and evaluation data. 
Therefore, it is essential to clarify the divergence between the 
two distribution trends for application. 

In the present work, the cross sections for the 
?7A1(n,2n)°°Al reaction in the 14 MeV neutron energy region 
are determined by utilizing the AMS technique with careful 
consideration of the experimental background to meet the re- 
quirements of nuclear physics and nuclear engineering. Af- 
ter comparing the (n,2n) cross sections of °7A1 obtained in 
the present measurement with the existing data and evalua- 
tion values, the new experimental results determined a new 
excitation function more accurately. 


I. EXPERIMENTAL PROCEDURE 
A. Neutron source 


In the neutron energy range from 13.5 to 15.0 MeV, irradi- 
ations were performed at the K-400 neutron generator of the 
China Academy of Engineering Physics (CAEP). The neu- 
tron yield of the K-400 neutron generator is approximately 
(3-5) 10!° n/s. The neutrons were produced by a T(d,n)“He 
reaction with a D+ beam current of 250 A. The thickness 
of the solid tritium-titanium (T-Ti) target used in the genera- 
tor was 2.5 mg/cm”, and the energy of the incident D+ beam 
was 300 keV. An Au-Si surface barrier detector was installed 
at 135° with respect to the D* beam and at a distance of 100 
cm from the T-Ti target for the detection of accompanying 
particles. The neutron flux was recorded and kept every ten 
seconds to obtain the neutron flux fluctuation over the entire 
irradiation period. 


B. Samples irradiation and neutron energy 


High purity aluminum (0.5 mm, 99.999% purity) foil was 
sandwiched between natural zirconium (0.1 mm, 99.5% pu- 
rity) and niobium (0.1 mm, 99.95% purity) foils. The first set 
of foils was tailored to form a round disk with a diameter of 20 
mm. Four sandwich samples were placed approximately 5.0 
cm away from the T-Ti target relative to the deuteron beam’s 
incident angles of 0°, 30°, 60° and 100°. The dimensions of 
the second set of foils were 0.5 cm in width and about 4 cm 
in length, which were in arc shape and placed about 3.0 cm 
away from the T-Ti target and about 60° to 135° relative to the 
deuterium ion incident direction. The schematic diagram of 
experimental geometry is shown in Fig. 2. The second set of 
foils was cut into 8 squares with a side length of 0.5 cm after 
the irradiation. The precise weight of each foil was deter- 
mined with an analytic balance (+0.01 mg) after irradiation. 

Two sets of aluminum samples placed at different locations 
from the tritium titanium target were irradiated for 20h. It is 
essential to accurately determine the neutron energy for cross- 
sectional measurements, particularly those near the thresh- 


Fig. 2. (Color online) Schematic diagram of experimental geometry. 


old. The neutron energies were measured using the Zr-Nb 
method [19, 20]. As a cross-check, the neutron energies of 
the first set of foils were calculated using an improved theo- 
retical method, which was used to determine the mean D-T 
neutron energies in a large solid angle sample [21, 22]. There 
are two sets of parameters for the Zb-Nb method, and the 
neutron energies were calculated using these two sets of pa- 
rameters. The mean value of the two results was used to de- 
termine the final neutron energy, as described in Sect. V. The 
uncertainties in the incident neutron energies given above are 
the quadratic roots of the summed square of the uncertain- 
ties caused by the results of the activity measurements for Zr 
and Nb foils [23]. The entire available energy range starts 
from the threshold of the (n,2n) reaction to approximately 
14.90 MeV. It is reasonable to suppose that the scattered neu- 
trons from the target construction and the surrounding en- 
vironment are insignificant because of the high threshold of 
the 27 A1(n,2n) 26A] reaction. In addition, the cross-section of 
°3Nb(n,2n)°2”Nb was used as a standard to monitor the inte- 
gral neutron flux for all samples. The relevant nuclear param- 
eters for °*Nb(n,2n)?~”Nb and ?°Zr(n,2n)*"”"*8Zr are listed in 
Table 1. 


C. Measurement of "Nb and *”"*®Zr activity 


After irradiation, there was a cooling period. Then, the 
activities of those Nb and Zr samples were determined us- 
ing an HPGe y-ray detector(as shown in Fig. 3). Before Nb 
and Zr samples were measured, a series of standard point 
sources (Na, Co, !33Ba, 137Cs, and 152pu) of known ac- 
tivity were used to determine the absolute full-energy peak 
efficiency of a lead-shielded high-purity germanium detec- 
tor (ORTEC, USA). The radioactivity of the reaction product 
was determined via the well-calibrated HPGe detector with 
a relative efficiency of 68% and an energy resolution of 1.82 
keV FWHM at 1.33 MeV. The distance between the sample 
and the detector was 5.0 cm. Efficiency calibration was de- 
termined using point-like standard gamma-ray sources. The 


detection efficiencies for the point source placed at a distance 
of 5.0 cm from the detector were determined using Eq. (1) 


C’ 


P~ Ao exp—*# AtL, (1) 
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where C” is the count number of the point source during the 
counting time (At), Ao is the source activity at the time of 
manufacturing, t is the time elapsed from the date of manu- 
facturing to the start time of counting, is the decay constant, 
and I, is the decay y intensity. 
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Fig. 3. (Color online) Typical y-ray spectrum in zirconium and nio- 
bium sample. 


The efficiency values obtained for the characteristic y-ray 
energies of the point source are described using an exponen- 
tial function [24] expressed in Eq. (2). The detector efficiency 
for the characteristic y-ray energies of °””"Nb and °°”"*8Zr was 
then obtained using Eq. (2). The detector efficiency uncer- 
tainties of °””"Nb and *””"*8Zr were obtained by evaluating the 
uncertainties of the standard sources. 


e(E) = eo exp(—E/Eo) + €e (2) 


The fitting parameter values are 2.53 (co), 264.14 (Eo), and 
0.33 (ec), respectively. The value of R? is 0.992. 


D. 7Al measurement at GYIG 1 MV AMS 


Because *°A has a half-life of millions of years, it is diffi- 
cult to measure 7°Al with high-purity germanium; therefore, 
an ultrasensitive AMS was needed to measure 7°Al in the 
present work. As a dedicated accelerator mass spectrome- 
ter, the GYIG 1 MV AMS instrument has been installed at 
the Institute of Geochemistry, Chinese Academy of Sciences 
since 2022 and was designed for multinuclide measurements 
such as Be, na. 26A, and “!Ca. After the installation was 
completed, the relevant acceptance and commissioning works 
were also completed. The present study was performed us- 
ing a 1 MV AMS. A schematic of the GYIG AMS facility 


TABLE 1. Associated decay data of °*Nb(n,2n)”"Nb, °°Zr(n,2n)"*8Zr and 77Al(n,2n)?°Al 


Reaction Abundance of target Q value Half-life of y-ray energy y-ray intensity 
isotope 7 (%) (MeV) product T; 2 (keV) I, (%) 

°3Nb(n,2n)””"Nb 100 —8.831 (10.15 + 0.02) d 934.44 99.15 + 0.04 

7 r(n,2n)"*8Zr 51.45 —11.965 (78.41 + 0.12) h 909.15 99.05 + 0.03 

21 A1(n,2n)” A1 100 —12.275 (7.14 + 0.24) x 10° a 1808.65 99.76 + 0.04 


(model 3SDH-1 Pelletron) is shown in Fig. 4 [25], which 
mainly includes a cesium sputter negative ion source, a low- 
energy injection system (electrostatic analyzer and LE mag- 
net), a tandem accelerator, a high-energy analysis system (HE 
analyzing magnet and electrostatic analyzer), and an energy- 
sensitive detector [26, 27]. 


AlO; was chosen as the proper sputter material. After 
neutron activation, approximately 0.08 g of aluminum foil 
was weighed and dissolved in 5 ml of hydrochloric acid (4.5 
M) to ensure a uniform mixture of the irradiated sample ma- 
terial. The solutions were diluted with deionized water to 20 
mL, respectively. The pH of the solution was adjusted by 
ammonium hydroxide to 9. A saturated ammonium oxalate 
solution was added to precipitate aluminum as aluminum hy- 
droxide (Al(OH)3). The Al(OH)3 precipitate was separated 
by centrifugation at 10000 rpm for 5 min, washed once with 
saturated ammonium oxalate solution, and rinsed twice with 
deionized water. The aluminum hydroxide precipitate was 
transferred to a ceramic crucible. The sample was dried at 
60 °C and then heated to 900 °C in a muffle furnace to ob- 
tain aluminum oxide powder. After mixing with Ag powder 
(Mass: Mass=1:1), the prepared samples were pressed into a 
copper holder, which acted as the cathode in the cesium sput- 
tering negative-ion source, capable of holding 40 cathodes of 
solid samples [28, 29]. 


Before the measurement of *°Al, the pilot °7Al~ beam was 
injected into the accelerator after being selected by the in- 
jection magnet, and the current of 7’Al~ was 0.1-1 A. Be- 
cause the electron affinity of aluminum is relatively low, the 
intensity of the °7A1~ current from the cesium sputtering neg- 
ative ion source is insufficient compared to C~% and BeO™ in 
this machine. Fortunately, Mg™ ions are unstable; thus, iso- 
baric interference is insignificant in the 26A] measurement. 
The accelerator terminal voltage was set to 0.41 MV, and 
negative ions were accelerated in the first part of the accel- 
erator. After passing through argon gas at a pressure of 27 
Torr, the negative ions were stripped into positive ions with 
different charge states via collisions and Coulomb breakups 
and then accelerated in the second part of the accelerator. 
°71'* particles were chosen by high-energy analysis system. 
The beam transmission was approximately 30%, which is the 
ratio of the 7’Al current measured in the Faraday cups be- 
fore and after passing through the high-energy magnet. After 
the 7’Al'* ions reached the Faraday cup before reaching the 
detector, all parameters were adjusted to fit the 26A] beam 
transport. Then, the sequential injection of 7’Al and *°A1 is 
achieved with a magnet bias supply coupled to the “bouncer 


magnet.’ The current of the °7a1'* ion was integrated at the 
offline Faraday cup, whereas the radioisotope 7°Al was mea- 
sured using an ionization chamber, which provides high res- 
olution in the energy measurement [30]. Under the influence 
of the sample preparation level, mixing uniformity with the 
conductive medium, and sample bombardment degree in the 
ion source, the multiple measurement current value of 27 A1" 
in a high-energy Faraday cup can fluctuate between tens and 
hundreds of nA, and the 26A] count follows the same pro- 
portion. Consequently, the %%A1/7 AI ratio measured between 
different rounds exhibited little deviation(as shown in Fig. 5). 
Because AMS is a relative measurement, the standard sample, 
processing the blank sample, and unknown sample should be 
measured sequentially. Meanwhile, the raw ratio results were 
subtracted from the processing blank value and normalized to 
the standard sample nominal value. 


To validate the performance of the newly installed instru- 
ment, a series of standard samples with different gradients 
and blank samples were used. The ratio of the blank sam- 
ples was less than 4 x 10715. The ratios of the Al stan- 
dard (5.00 x 1071!) were reproduced to be better than 1% 
(as shown in Fig. 6). To measure the (n,2n) cross-section, 
the raw ratios(before normalization) of 26 A1/7 Al are between 
6 x 107! and 10718, which are much larger than the ma- 
chine background of 7°A1/*7Al at the GYIG 1 MV AMS. To 
verify the measurement results of the GYIG 1 MV AMS, four 
samples and one blank sample were sent to the Xi’an Accel- 
erator Mass Spectrometry Center. The measurement standard 
for AI” Al is 1.065 x 107"! in Xi’an and the background 
is 4.71x1071! [30]. The values of the four irradiated sam- 
ples were between 1.771071! and 2.09 10~ 1", which is in 
good agreement with the results of GYIG 1 MV AMS. 


Il. DATA ANALYSIS 


The experimental cross-section Gexpt for the 21 A1(n,2n) 6A] 
reaction can be determined using Eq. (3) [27]: 


N26 1 
Noz Diotal 


(3) 


Oexpt = 


The ratio %A1/7 A1 was determined by AMS, as shown in 
Table 3, and otai is the neutron fluence that can be measured 
using a monitor. 


Fig. 4. (Color online) Schematic of the GYIG 1 MV AMS, Model 3SDH-1. I: Ion source II: Sequencer and LE mass analyzer. HI: 1 
MV Pelletron accelerator and gas stripper. IV: HE analyzing magnet, off-axis Faraday cups. V: The rare isotope beamline (containing the 
retractable stopper foil, the electrostatic spherical analyzer, the 45° analyzing magnet). VI: Gas ionization detector [25]. 
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Fig. 5. (Color online) *°Al counting,” Al'* current and 7°Al/?’Al 
ratio fluctuation among different rounds. 


A. Neutron fluence 


The neutron fluence of the aluminum sample during irradi- 
ation was calculated using Eq. (4) [23]. 


FOXA 


MNanlyeKSDom ’ (4) 


Protal = oT = 


where ¢ is the neutron flux averaged over the sample in the to- 
tal irradiation time 7’, cm is the monitor reaction cross section 
for °°Nb(n,2n)°”"Nb, S = 1—exp~*7 is the growth factor of 
the residual nuclide, A is the decay constant of 92"NIb, cis the 
full-energy peak efficiency of the 934.44 keV y ray, I, is the 
emission probability of the 934.44 keV y ray in "Nb decay, 
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Fig. 6. (Color online) Raw measured isotope ratios for various ref- 
erence samples versus their nominal values. 


n is the isotopic abundance of the target nuclide Nb, M is 
the mass of Nb, D = exp~*“! — exp~”2 is the counting 
collection factor, tı (t2) is the time intervals between the end 
of irradiation and the start (stop) of counting, A is the atomic 
weight of Nb, C is the measured full-energy peak area, NA 
is the Avogadro constant, F is the composite correction factor 
of the activity expressed by Eq. (5), and K is the neutron flux 
fluctuation factor expressed by Eq. (6) [31-33]. 


F = FF Fr, (5) 


where Fs, Fẹ and F, are the correction factors for the self- 
absorption of y-rays with a given y energy, counting geom- 


etry, and coincidence summing effect of the cascade y-rays 
in the investigated nuclide, respectively. In this study, all the 
samples were measured 50 mm from the surface of the detec- 
tor. This distance is so large that F, and F, are considered 
negligible [32]. The self-absorption of 934.44 keV y-ray by 
0.1 mm Nb was calculated to be 1.003. 


Ka Li bill = ieee exp 47 © 


where L is the number of time intervals, in each of which the 
accompanying œ counts were individually recorded, t; is the 
duration of the ¿™ time interval, T; is the time interval from 
the end of the 7“ interval to the end of irradiation, Q; is the 
relative neutron flux averaged over the sample during t; (i.e., 
a counts in t;), ġ is the relative neutron flux averaged over 
the sample in the total irradiation time T (i.e., total œ counts 
in T) [34]. 


B. Uncertainty analysis 


The total uncertainty of the °7A1(n,2n)*°Al reaction con- 
sists of the monitoring reaction and the AMS contributions. 
Individual contributions to the total uncertainty are listed in 
Table 2. The uncertainty in the AMS measurements includes 
the statistical uncertainty (4.9%-50.2%), reproducibility of 
the measurement (1.0% for the standard sample), and sys- 
tematic contribution from the measurement relative to refer- 
ence materials (1.0%). As the most commonly used standard 
cross-section, the uncertainty of °*Nb(n,2n)°~”Nb is approx- 
imately 1.5% [5]. The uncertainty (3.0%) of the HPGe effi- 
ciency calibration was derived from the accuracy of the cal- 
ibration sources and calibration procedure. The total uncer- 
tainty was estimated using the square root of the quadratic 
sum of the above sources, and the overall uncertainties ranged 
from 7% to 50.2%. 


IV. THEORETICAL CALCULATIONS 


To obtain accurate and trustworthy nuclear data, theory, 
and experiments must be mutually verified [34-36]. Based on 
this, theoretical calculations for the 21 A1(n,2n)*°Al reaction 
were performed using the statistical nuclear reaction model 
codes EMPIRE-3.2.3 and TALYS-1.95 [37]. Many studies 
have focused on introducing these two codes, both at home 
and abroad [3 1-36], which will not be repeated in this study. 
Here, a brief introduction of the basic principles is provided. 
Computational nuclear reaction mechanisms vary depending 
on the energy of the incident particle [35]. These algorithms 
include several nuclear models that use distinct sets of optical 
model parameters and level densities to determine the con- 
tributions of the direct reaction (DI), pre-equilibrium emis- 
sion (PE), and compound nucleus (CN). To perform theoret- 
ical calculations, the best possible input parameters must be 
used, which are considered to replicate the most satisfactory 
outcomes compared to the currently obtained data and all pre- 
vious experimental data [38, 39]. In this study, cross-sections 


for the 21 A1(n,2n) 26A] reaction were calculated using the de- 
fault parameters with an optical model (EMPIRE-3.2.3) and 
different level density models (TALYS-1.95). The calculated 
outcomes and evaluated databases are shown in Fig. 4. 


V. RESULTS AND DISCUSSION 


The results of the neutron energy and cross sections of the 
27 A1(n,2n) 26A] reaction in the 14 MeV region are summarized 
in Table 3 and Fig. 7. In Fig. 7, the cross sections of the 
°7A1(n,2n)*°Al reaction are compared with the available lit- 
erature data in the EXFOR database. The measurements of 
13-15 MeV neutron energies showed discrepancies between 
the data from different groups. Figure 7 shows a good agree- 
ment with the data of Wallner et al.[5] and Filatenkov et al., 
obtained within the experimental uncertainties. In addition, 
the experimental results of Zhao et al. [14] are inconsistent 
with the results of this study. Below the neutron energy of 
14.0 MeV, the data of Sasao et al. [12] are consistent with the 
present data. The data of Iwasaki et al. are two times higher 
than the present work, between 14.0 MeV and 15.0 MeV. 
Compared to gamma spectrometry, the AMS technique of- 
fers a more competent method for the measurement of cross- 
sections of reactions with long-half-life product nuclides. Be- 
cause of the short irradiation time, the statistical uncertainty 
of 7°Al became the primary source of the uncertainty source 
term in the experiment. If sufficient irradiation time is pro- 
vided, the statistical uncertainty is greatly reduced. Without 
considering the AMS statistical uncertainty, the accuracy of 
this experiment was significantly improved. Therefore, the 
present results clarify the divergence of the two trends among 
previous experimental data, and a conclusion can be drawn 
as follows: the results that have a consistent trend with the 
present data are more accurate. 
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Fig. 7. (Color online) Comparison of reported cross sections of 
?7Al(n,2n)°Al reaction from threshold to 20 MeV region. 


The evaluation values compiled by BROND-3.1, JEFF- 
3.3, ENDF/B-VIII.0, CENDL-3.2, and JENDL-5.0, are also 
shown in Fig. 7. The measured 27 Al(n,2n)°Al reaction cross 


TABLE 2. The uncertainty source analysis of measured 7’Al(n,2n)”°Al reaction cross section 


Source of uncertainty Uncertainty (%) 
Counting statistics of *°Al 4.9-50.0 
Reproducibility of standard sample ~1.0 
Counting statistics of °””"Nb 1.0 
Half-life of °”"Nb 0.2 
Sample mass 0.1 
Neutron flux fluctuation factor [19, 32] 3.0 
Cross section of °*Nb(n,2n)?2”Nb [5] ~ 1.5 
Standard sample (7°AI/7’Al) 1.0 
Total correction factor (F) 1.0 
Branching ratio 0.5 
Detection efficiency 3.0 
Total 6.9-50.2 


TABLE 3. The present cross-section results of the 7’Al(n,2n)”°Al reaction at 14 MeV neutron energy region 


Measured Neutron energy Neutron fluence The ratio of Cross section 
ratio (MeV) (n/cm?) SAI” AI (mb) 
1.85 + 0.05 14.91 + 0.18* 5.62 x 10"? (2.09 + 0.10) x 10713]* 37.20 + 2.80 
1.75 + 0.05 14.80 + 0.16* 3.50 x 10"? (1.07 0.07) x 10713]¢ 30.60 + 2.72 
1.60 + 0.05 14.51 + 0.15* 3.77 x 10” (7.18 + 0.59) x 10714]* 19.00 + 1.92 
1.54 + 0.05 14.37 + 0.09 1.22 x 1018 (1.69 + 0.14) x 10713]? 13.79 + 1.42 
1.40 + 0.05 14.17 + 0.09 1.44 x 1018 (1.46 + 0.13) x 10713]? 10.13 + 1.10 
1.35 + 0.05 14.10 + 0.09 1.39 x 1078 (1.01 + 0.16) x 10713]? 7.27 + 1.20 
1.31 + 0.04 13.99 + 0.09* 3.46 x 10’? (1.77 + 0.32) x 10714]¢ 5.11 + 0.97 
1.29 + 0.04 13.95 + 0.09 1.46 x 1013 (4.52 + 0.89) x 10714]? 3.10 + 0.63 
1.23 + 0.04 13.86 + 0.08 7.10 x 10”? (1.60 + 0.38) x 1074]? 2.26 + 0.55 
1.11 + 0.04 13.70 + 0.07 9.68 x 10"? (7.72 + 2.24) x 10719]? 0.80 + 0.24 
1.09 + 0.04 13.68 + 0.07 1.25 x 1013 (4.47 + 1.83) x 1077]? 0.36 + 0.15 
1.03 + 0.04 13.61 + 0.06 1.42 x 1018 (4.09 + 2.05) x 10719]? 0.29 + 0.15 


Note: In the Second column, * means that the corresponding neutron energy was calculated by the theoretical method; In the third column, 


superscript a means they are the first set of samples, and b means that they are the second set of samples. 


section falls within the range of values calculated by differ- 
ent models at 14 MeV neutron energy range. Our values 
were almost in agreement with those of CENDL-3.2. The 
near-threshold region above 13.5 MeV is based on different 
statistical model calculations and thus differs from evalua- 
tion to evaluation. The datasets of TENDL-2021 are identical 
and use the recent nuclear model code TALYS-1.95, where 
the parameters have been adjusted to reproduce the exist- 
ing experimental data; however, because there is not enough 
experimental data to support the evaluated cross-section in 
the 15-20 MeV neutron energy region, there is still a di- 
vergence among the energy-dependent neutron cross-sections 


from the evaluated cross-section libraries JEFF-3.3, CENDL- 
3.2, ENDF/B-VIII.0, and JENDL-5.0. Therefore, further ex- 
perimental studies are required. 


In addition, the excitation function of the 21 A1(n,2n)*°Al 
reaction was calculated theoretically using TALYS-1.95 code 
and EMPIRE-3.2.3, with default parameters. In Fig. 8, the 
shapes of the excitation curves calculated using TALYS-1.95 
(Idmodels 1-6) and EMPIRE-3.2.3 exhibit a trend similar to 
the present data set, which increases with increasing neutron 
energy around 14 MeV. However, the results calculated using 
Idmodel 1-6 in TALYS-1.95 significantly overestimated the 
excitation function. The excitation functions calculated by 
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Fig. 8. (Color online) Excitation function for the 21 Al(n,2n)™Al re- 
action from threshold to 20 MeV region. 


the level density model Idmodel 1-2 in TALYS-1.95, are in 
good agreement with the experimental data compared to the 
level density model Idmodel 3-6. From the threshold to 15.5 
MeV, EMPIRE-3.2.3 is consistent with the results reported by 
TENDL-2021 within the uncertainties. Therefore, more high- 
precision experimental data are still needed to provide param- 
eters and an experimental basis for theoretical model calcula- 
tions, especially for energy above 15.5 MeV. The present re- 
sults contribute to improving the knowledge of cross-sections 
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